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Abstract 

We consider corrections to bino annihilation due to sfermion mixing. 



Within the standard model of electroweak interactions, the only plausible dark matter 
candidate is the neutrino. A neutrino (v T , for example) with a mass of 10-30 eV would make a 
significant contribution to the overall mass density of the Universe. A light neutrino, however, 
is a hot dark matter candidate and as such cannot account for the observed structure on small 
(galactic) scales. Cold dark matter alone also has its difficulties. It is difficult to make large 
scale structure and conflicts with observation of the microwave background anisotropies. At 
present, the best mix of dark matter and baryons seems to be 60-70% cold, about 30% hot, 
and the remainder in baryons (see |l| and references therein) . Cold dark matter still appears 
to be a necessity. Before their demise, heavy (masses >, a few GeV) neutrinos were by 
far the simplest cold dark matter candidates. The minimal supersymmetric standard model 
(MSSM) is perhaps the simplest logical extension which offers a cold dark matter candidate. 
Because of an unbroken discrete symmetry, R-parity, the lightest supersymmetric partner 
(LSP) is expected to be stable. This makes the LSP a natural candidate for role of the 
non-baryonic dark matter. 

The identity of the LSP is concealed by the relatively large number of unknown param- 
eters of the MSSM. Nevertheless, by a combination of i) theoretical presumptions (eg that 
the MSSM will eventually be contained in a grand unified theory), ii) limits from acceler- 
ator experiments and iii) cosmological bounds, it is possible to identify the more plausible 
regions of the parameter space. In previous studies, the lightest neutralino, a linear com- 
bination of the supersymmetric partners of the neutral SU(2) gauge boson (neutral wino), 
£7(1) hypercharge gauge boson (bino), and the two neutral components of the Higgs doublets 
(higgsinos) has been suggested as the most likely candidate for the LSP 0]. The mass matrix 
for these states contains three unknown parameters: (1) M 2 , the supersymmetry breaking 
SU{2)l gaugino mass; (2) e, the higgsino mixing mass; and (3) tan/5 = t>i/t>2, the ratio of 
the vacuum expectation values of the Higgs doublets. In a large portion of this parameter 
space, the least massive eigenstate is a nearly pure bino ||. 

Restricting our attention to the bino B as the candidate LSP, it is straightforward to use 
its couplings to the other fields in the MSSM to compute the relic abundance. These cal- 
culations were previously performed under a variety of assumptions: degenerate squark and 



1 



slepton masses with m~ = mj > mg taken as a relatively free parameter [|, [5] ; unequal 
squark and slepton masses with values run down from the GUT scale using renormalization 
group equations with supergravity-inspired boundary conditions, and sfermion mixing ne- 
glected H ; the same with sfermion mixing included [|7|, |§ . Note there is also a great variety 
in the level of detail of these calculations. Only in |5], || have all LSP-LSP annihilation 
channels been included. In || the annihilations of the LSP and the next lightest neutralino 
(co-annihilations |§) were also included. 

In general, the sfermion mass matrix contains off diagonal terms connecting the left and 
right sfermions. These are proportional to the fermion-Higgs coupling and hence to the 
fermion mass. Sfermion mixing is also dependent on the higgsino mixing mass e. Since the 
region of the parameter space containing the almost pure bino is in the range of large e, for 
mg > m t we expect the off-diagonal terms to have a large effect on bino annihilation by 
exchange of a top squark. 

Though mixing was included in ||, all parameters were set according to a particular 
GUT renormalization (albeit a good one). The purpose of this paper is to examine in more 
detail the full effect of including sfermion mixing. A subsequent paper will deal with another 
correction to the bino annihilation process, s-channel annihilation to fermions via a Z° at 
the one loop level. 

We begin by deriving the annihilation cross-section for binos, keeping the explicit depen- 
dence on sfermion mixing. The general form of the sfermion mass matrix is [H| 

_ _ (M 2 L m 2 \ (f L \ 
Ul (J Mi) (/J ' (1) 

where m 2 = rrif(Af — ecot/?) for weak isospin +1/2 fermions and rrif(Af — etan/3) for 
weak isospin —1/2. Af is the soft supersymmetry breaking trilinear mass term arising from 
superpotential terms such as Hi y 2Q{u c ,d c ) or H2Le c , and e is the Higgs mixing parameter, 
often written as — //; see [|11|] for an overview of couplings in the MSSM. This mass matrix is 
easily diagonalized by writing the diagonal sfermion eigenstates as 

fi = hcosOf + f R sm9 f , 

h = -f L sm9 f + f R cos9 f . (2) 
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With these conventions we have the diagonalizing angle and mass eigenvalues 

|-| arctan |2m 2 /(M 2 - M|) 
±arctan|2m 2 /(M| - M|)|, 



7/ = sign[— m 



M 2 >M R , 



M 2 <M 2 R , 



rn 



1,2 



- (M 2 R + M 2 L ) T J(M 2 R -M 2 y + 4m* 



(3) 



Here is chosen so that mi is always lighter that m 2 . Note that in the special case 
M L = M R = M, we have 8f = sign[— m 2 ](7r/4) and m\ 2 = M 2 =f \m 2 \. 

From here it is straightforward (if tedious) to enumerate the Feynman diagrams contribut- 
ing to the annihilation process. The dominant contribution is due to sfermion exchange and 
is derived from bino-fermion-sfermion interaction lagrangian, 



C — 

ffB 



-^9' (Y R fP L Bf R + YjP R Bf L ) + h.c. 

±g'f (J x xP L + f 2 wP L + hyP R - f 2 zP R ) B + h.c. 



(4) 



where x = Y R sin8 f , y = Y L cos8 f , w = Y R cos8 f , z = Y L sm8 f , P R)L = (l±7 5 )/2, Y R = 2Q f 
and Yl = 2(Qf — T 3 f) where Qf is the fermion charge and T 3 j is the fermion weak isospin. 

To derive a thermally averaged cross section, we make use of the technique of ref. |T2 
We expand (crv re i) in a Taylor expansion in powers of x = T/m^ 



\CVrei) = a + bx + 0(x 2 



The coefficients a and b are given by 



-3 + 



3m 2 , 

Av 2 f m 2 ~ 
J b , 



a .f 



(5) 



(6) 



(7) 



where df and 6/ are computed from the expansion of the matrix element squared in powers 



of p, the incoming bino momentum, and Vf 
space integrals. 



m 2 f/m 



B' 



is a factor from the phase 
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We summarize the result by quoting the computed a/: 
~ 3" \ 

; 1287T 

where Aj = m 2 ? + m~ — mj. The result for bf is too lengthy for presentation here, but was 
computed and used in the numerical integrations to obtain the relic densities. The results 
reduce to the results quoted in |3j in the limit of no mixing and equal left and right sfermion 
masses. 

We next discuss the effect of considering this mixing interaction on the relic density and 
detectability of binos. In all of the following cases we have chosen a value of tan (3 = 2 and 
e = ±3000 GeV and assume a top quark mass of 120 GeV. This choice of e places us in the 
region of parameter space where the LSP is a nearly pure bino for a wide range of gaugino 
masses, M2 — 50 — 5000 GeV. We also take the diagonal part of the sfermion mass matrix 
to have a constant value for all squarks and sleptons, and denote the common value of Ml 
and Mr by M. This assumption works best when the sfermion masses are large compared 
with m t , nig, and Mz- For simplicity, we take this condition to apply in all cases. Previous 
results can be viewed as a special case of this more general form. Indeed, when mixing is 
neglected one must make a particular choice of Af = ecot (3 and Aj = etan/3 for T 3 = ±1/2 
fermions. Thus results for the relic bino density when mixing is neglected implicitly assume 
that a specific (and different) value of Af is chosen for each value of tan/? and e. In the 
mixed case we will also take a common value for all Af and denote it by A. 

Sfermion mixing has several effects on the cross-section that can be readily understood. 
First, the cross-section is no longer p-wave suppressed. In (||), we see that the zero temper- 
ature cross-section now has a piece proportional to mg 2 /A 2 rather than m 2 /A 2 . (Without 
mixing, when ni\ = 777.2, this piece disappears.) For mg > rrif, this serves to increase the 
cross-section. However this aspect of mixing produces little effect on b, which already in- 
cludes terms proportional to mg 2 /A 2 . On the other hand, the splitting of the sfermion 
masses makes A2 > Ai. The cross-section (relative to the unmixed case) may be either 
larger or smaller depending on which mass term is held fixed. When the diagonal elements 
are fixed, the presence of the lighter sfermion in the mixed case tends to increase the cross- 
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Ai(m/w — 2rrigWZ + rrifZ ) ± A2(m/x ± 2m^xy + rrify z 

AiA 2 



section; however, since the bino couples with different strengths to the two sfermion mass 
eigenstates, the cross-section can actually decrease [see fig. (1)]. When the lightest mass 
eigenvalue mi is fixed, the cross-section goes down as half the sfermions get heavier. 

In fig. (1), we show the effect of mixing on the coefficient a given by @ as a function 
of the bino mass (which for our chosen values of e is directly proportional to the gaugino 
mass M2). In the upper and lower curve where mixing is included, we have taken Af = 
and chosen M so that mi = mg, the smallest mass consistent with the identification of the 
bino as the LSP. In the upper curve e = 3000 GeV, while in the lower curve e = —3000 GeV. 
The middle curve, where mixing is neglected, has the same value of M as the other curves 
but with off diagonal mass entries set to zero by adjusting the Af's to the appropriate 
values. In these extreme cases, mixing may either increase or decrease the cross-section 
depending on the sign of e. The increase in cross-section is primarily due to the presence of 
a lighter sfermion in the propagator. When e = —3000 GeV, there is a cancellation which 
then decreases the cross-section at larger mg. (However, in this case the relic abundance 
turns out to be too high, as we will discuss below.) Notice also the large effect of crossing 
the top quark threshold. 

In fig. (2), we fix the relic density Q^h 2 and consider the effects of mixing on the zero 
temperature cross-section. fi~ is the mass density of the binos in units of the critical density, 
and h is the Hubble parameter in units of 100 km s -1 Mpc -1 . In the mixed case, we set 
Af = 0, and for each value of mg, vary M until the value of fl^h 2 equals 1/4, corresponding 
to a critical mass density of binos with a value of h = 1/2. Here e = 3000 GeV only, as there 
are no solutions to Q^h 2 = 1/4 with e = —3000 GeV. In the unmixed case, we set the Af's 
appropriately and again adjust M so that Q^h 2 = 1/4. The ratio of the a values in the 
mixed and unmixed cases is plotted in fig. (2). Notice here that a is approximately twice as 
large in the mixed case. 

The enhanced value for the zero-temperature cross-section has implications for the pros- 
pects of detecting neutralino dark matter. Among the proposed effects signaling the presence 
of the particles comprising dark matter are (i) gamma ray line flux from annihilations in the 
galactic halo [[H| (ii) high energy neutrino flux from dark matter particles annihilating after 
being trapped in the Sun or the Earth [111] and (iii) direct cryogenic detection of DM particles 
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scattering off nuclei []I3J (for a review, see Referring again to fig. (2) we see that the 

inclusion of mixing allows the relic abundance of binos to be the same as in the nonmixing 
case in spite of the increased cross section a, due to a large decrease in the p-wave part of the 
cross-section given by b. The increased value of a will directly increase the expected gamma 
flux in (i) as the flux is proportional to a. Bino annihilation in the sun is not affected, as 
the rate there is limited by the sun's ability to trap binos. In cryogenic detectors, the effect 
will also be important, as pointed out in [lTj, due to the effect of mixing on the elastic 
cross-section (see also [HJ for more on the elastic scattering of binos). 

In fig. (3), we show the region of the mg — A parameter space excluded by the conditions 
1) Qgh 2 < 1/4, 2) the lightest sfermion is heavier than the bino, and 3) the lightest sfermion 
is heavier than 74 GeV. Here we have set e = 3000 GeV, but note that for other choices of 
e the only substantial change in the picture is a shift in the values of Af (since the stop 
is dominant over much of the space, the relic abundance depends largely only on the off- 
diagonal matrix elements for the stops). The upper limit on the bino mass has the same 
physical origin as in the unmixed case. For no mixing, the smallest value of occurs when 
= mj ; as both mg and mj are increased, f2g also increases. When we include mixing, 
the upper limit depends on the trilinear mass term A, as shown in fig. (3). Again, the 
smallest value of for a fixed Af and mg occurs when the bino is as heavy as the lightest 
sfermion. 

The lower limit on mg seems to be a general feature of including mixing. The stop mixing 
effect becomes large outside of a region a few hundred GeV wide near Af — 1500 GeV and the 
condition that mg > mj forces a large value for the diagonal entries, M. The large sfermion 
masses suppress the annihilation rate below the top threshhold. Near A = 1500 GeV, there 
is little stop mixing, M can be much smaller, and binos can efficiently annihilate into r's and 
6's. In this region, the condition mj > 74 GeV limits the annihilation rate and determines 
the lower bound on mg. 

Finally, in fig. (4) we show the effect of A in the off-diagonal mass term. The bino mass 
is fixed at 150 GeV, e = 3000 GeV and we plot Q^h 2 as a function of A for several values of 
M. We effectively pass through the nonmixing case at a value of A = 1500 GeV (where the 
stops become unmixed). Again the points plotted are the largest range of A consistent with 
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the identification of the bino as the LSP. In the absence of mixing, as the sfermion mass 
was increased the cross section decreased leading to increased fl^h 2 . So choosing a desired 
value for fl^h 2 fixed the value of the sfermion mass in relation to the bino mass. Now with 
the addition of the parameter A it is possible to achieve a range of possible sfermion masses 
consistent with a given relic density as is evidenced in fig. 4. 

To summarize, we have calculated the effect of right and left sfermion mixing on the 
annihilation cross-section for binos. For a fixed galactic halo density of binos, we find mix- 
ing enhances the zero-temperature cross-section; this may improve prospects of detecting 
neutralino dark matter through gamma ray flux from annihilations in the halo or cryogenic 
detection of binos scattering from nuclei. When the mixing effect for stops becomes pro- 
nounced we also find stronger upper and lower bounds on mg coming from the requirement 
that Qgh 2 < 1/4 and the consistency condition that the bino be the LSP. Finally, given a 
halo density fig/i 2 , bino mass mg and choice of e we note that in the mixing case we are 
allowed a range of sfermion masses through variation of the trilinear mass terms Af, whereas 
the sfermion mass was uniquely determined in the unmixed case. 
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Figure Captions 



Figure 1: The zero-temperature cross-section, a, as a function of the bino mass when 
the diagonal entries of the sfermion mass matrix are held fixed. The upper 
curve (e = 3000 GeV) and lower curve (e = —3000 GeV) assume sfermion 
mixing. In the middle dotted curve the values of Af are taken in such a way 
as to create no mixing, but with the same diagonal mass entries as in the 
mixed cases. 

Figure 2: The ratio of the zero-temperature cross-sections, a m i X /a nom i X as a function of 
the bino mass when the relic density Vt^h 2 = \ is fixed, e = 3000 GeV and 
the diagonal mass choices are made as in Figure 1. 

Figure 3: Allowed regions (Qgh 2 < j) in the mg — Af plane for e = 3000 GeV. 

Figure 4: The relic bino density il^h 2 at e = 3000 GeV as a function of Af for several 
values (as labeled) of the diagonal entry of the sfermion mass matrix, M. 



10 



10" f < 1 < 1 < 1 < 1 < r 




100 200 300 400 500 600 



m (GeV) 



This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9308324vl 



1.9 




-| g I I I I I I I I I I 

1 20 140 160 180 200 220 



m (GeV) 



This figure "figl-2.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9308324vl 




A (GeV) 



This figure "figl-3.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9308324vl 




A (GeV) 



